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ABSTRACT 

Context. Ultraviolet and optical surveys are essential to gain insight into the processes driving galaxy formation and evolution. The 
rest-frame ultraviolet emission, redshifted to optical and near-infrared for high-redshift galaxies, is key to measure the cosmic star 
formation rate. However, the ultraviolet light is strongly reddened and absorbed by dust. In starburst galaxies, the ultraviolet colour 
and the attenuation are intrinsically linked, allowing to correct for dust extinction. Unfortunately, evidence has been accumulating 
that the relation between ultraviolet colour and attenuation is different for normal star-forming galaxies when compared to starburst 
galaxies. 

Aims. It is still not understood why star-forming galaxies deviate from the ultraviolet colour-attenuation relation of starburst galaxies. 
Previous work and models hint that the role of the shape of the attenuation curve and the age of stellar populations have an important 
role. In this paper we aim at understanding the fundamental reasons to explain this deviation. 

Methods. We have used the CIGALE spectral energy distribution fitting code to model the far ultraviolet to the far infrared emission 
of a set of 7 reasonably face-on spiral galaxies from the Herschel Reference Survey on a pixel-by-pixel basis. We have explored the 
influence of a wide range of physical parameters to quantify their influence and impact on the accurate determination of the attenuation 
from the ultraviolet colour, and why normal galaxies do not follow the same relation as starburst galaxies. 

Results. We have found that the deviation from the starburst relation can be best explained by intrinsic ultraviolet colour differences 
between different regions in galaxies. Variations in the shape of the attenuation curve can also play a secondary role. Standard age 
estimators of the stellar populations such as the i)4000 index or the birthrate parameter prove to be poor predictors of the intrinsic 
ultraviolet colour These results are also retrieved on a sample of 58 spiral galaxies drawn from the Herschel Reference Survey sample 
when considering their integrated fluxes. 

Conclusions. When correcting the emission of normal star-forming galaxies for the attenuation, it is crucial to take into account 
possible variations in the intrinsic ultraviolet colour of the stellar populations as well as variations of the shape of the attenuation 
curve. 

Key words, galaxies: star formation, galaxies: spiral, ultraviolet: galaxies, infrared: galaxies 



1. Introduction 

Over the last decade, large, deep surveys of high redshift galax- 
ies have been dedicated to gaining insight into the physical pro- 
cesses at play in the formation and evolution of galaxies across 
the Universe. The major process in the transformation of bary- 
onic matter is star formation. It converts the local gas reservoir 
into heavy elements that are ejected from star-forming regions 
by way of feedback, seeding the interstellar medium and the 
intergalactic medium with metals. One of the main constraints 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



on cosmological models is the so-called cosmic star formation 
rate (SFR) density, which has been intensely studied since the 
seminal work of Madau et al. (1996) both observationally (e.g. 
Steidel et al. 1999; Hopkins 2004; Perez-Gonzalez et al. 2005; 
Schiminovich et al. 2005; Bouwens et al. 2009; Magnelli et al. 
2009; Reddy & Steidel 2009; Rodighiero et al. 2010; van der 
Burg et al. 2010; Magnelli et al. 2011) and theoretically (e.g. 
Kitzbichler & White 2007; Dave et al. 201 1). 

To measure the SFR across the Universe, the ultraviolet 
(UV) is theoretically the wavelength domain of choice for high- 
redshift galaxies. Indeed, rest-frame UV radiation is redshifted 
into optical and near-infrared bands which are easily accessi- 
ble with broadband observations from the ground. UV is a direct 
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tracer of star formation as it is sensitive to the photospheric emis- 
sion of massive stars. Unfortunately the presence of dust affects 
its effectiveness: it reddens the UV-optical spectral energy distri- 
bution (SED) as it absorbs energetic radiation that is re-emitted 
in the mid- to far-infrared (IR), and it is also affected by scatter- 
ing out and into the line of sight. Following the usual convention, 
extinction encompasses absorption and scattering out of the line 
of sight while attenuation also takes into account the scattering 
into the line of sight. Correcting the UV emission for the atten- 
uation is therefore a crucial requirement to estimate the actual 
SFR of galaxies. 

A powerful way to correct for the attenuation is to combine 
attenuation sensitive star formation tracers with the IR emission 
(e.g. Calzetti et al. 2007; Leroy et al. 2008; Kennicutt et al. 
2009; Hao et al. 2011). Unfortunately IR data are not neces- 
sarily available and often still not at a sufficient depth and/or 
resolution, even with Herschel. In this context, one of the most 
commonly used methods to correct for the attenuation is the so- 
called IRX-yS relation (Calzetti et al. 1994; Meurer et al. 1999; 
Calzetti et al. 2000). This method links the observed UV slope p 
to IRX, a measure of the attenuation: 

IRX^XogiLdustlLFuv), (1) 

with Ldust the total luminosity of the dust, and Lpuv the far UV 
luminosity computed as vLy {X) with v the frequency and Ly {X) 
the monochromatic luminosity per unit frequency at the wave- 
length A. While this relation provides us with accurate estimates 
of the attenuation for starburst galaxies (Meurer et al. 1999), 
it fails, sometimes considerably, for galaxies forming stars at a 
lower rate. The fundamental reason for this deviation from the 
starburst law is still poorly understood. Various studies hint at 
the possible role of the age of the stellar populations and/or vari- 
ations of the shape of the attenuation curve (Bell 2002; Kong 
et al. 2004; Burgarella et al. 2005; Seibert et al. 2005; Cortese 
et al. 2006; Gil de Paz et al. 2007; Johnson et al. 2007; Panuzzo 
et al. 2007; Cortese et al. 2008; Boquien et al. 2009; Munoz- 
Mateos et al. 2009; Wijesinghe et al. 2011). This is a severe 
problem as surveys go ever deeper and probe the fainter end of 
the luminosity function which is populated by galaxies that can 
strongly deviate from the starburst law. The now routine detec- 
tion of normal star-forming galaxies at high redshift shows that 
there is urgency to understand what is the physical origin of this 
deviation. 

Observations over a broad wavelength range is key to under- 
standing the physical processes at play in star-forming galaxies. 
Indeed, the energy output of young stellar populations is domi- 
nated by the photospheric emission of short-lived, blue, massive 
stars that emit in the UV. These stars also ionise the surrounding 
gas which recombines emitting lines that can dominate the flux 
in optical bands (Anders & Fritze-v. Alvensleben 2003). At the 
same time the optical is also a combination of young and old stel- 
lar populations. The relative weight of the emission from these 
two populations changes with wavelength with younger popula- 
tions dominating at short wavelengths whereas the evolved pop- 
ulation dominating at longer wavelengths, peaking in the near- 
IR. To disentangle these populations a full sampling in the UV, 
optical, and in the near-IR is needed. These bands are how- 
ever sensitive to the dust which absorbs the emission of stars 
and re-emits the energy at longer wavelength in the mid- and 
far-IR. The emission in the mid-IR arises principally from hot 
dust stochastically heated and emission bands from large com- 
plex molecules such as PAH. Conversely the far-IR is dominated 
mainly by warm dust (~ 50 K, under -100 yum) and cold dust 



(~ 20 K, beyond ~ 100 yum), as can be seen in M33 for in- 
stance (Kramer et al. 2010). The emission of the warm dust in 
the far-infrared is particularly important to constrain the frac- 
tion of extinguished star formation as it dominates the energy 
budget. Understanding why star-forming galaxies deviate from 
the starburst relation therefore requires a multi-wavelength data 
set from the UV to the FIR. 

In this article we investigate, disentangle, and quantify the 
role of a large set of physical parameters to explain the devia- 
tion of star-forming galaxies from the starburst IRX-yS relation. 
To do so, we study the IRX-yS relation using FUV to FIR data 
for seven resolved galaxies from the Herschel Reference Survey 
(HRS, Boselli et al. 2010). 

In Sec. 2 we explain how we selected the galaxies to be stud- 
ied and how we processed the multi-wavelength data. Each data 
point has been modelled using CIGALE (Code Investigating 
GALaxy Emission, Noll et al. 2009a) to reproduce their SED 
from the FUV to the FIR as presented in Sec. 3. In Sec. 4 we 
analyse the direct estimate of a large number of physical pa- 
rameters allowing us to perform a detailed study of the relation 
between these parameters, and IRX, and in particular to con- 
strain the effect of the age of the stellar population and the shape 
of the attenuation curves. Finally, we conclude in Sec. 5 

2. Sample and data 

2.1. Sample selection 

The sample selection is guided by the core question of this paper: 
why and how do normal spiral galaxies forming at most a few 
Mq yr"^ of stars deviate from the IRX-/3 relation that has been 
derived for starburst galaxies? To answer this question we must 
be able to disentangle in a resolved way the emission from the 
different components of the galaxy: the young and evolved stel- 
lar populations, the nebular emission, and the dust. This yields 
constraints on the required bands and on the resolution. In addi- 
tion to ensure that this study is not affected by systematic vari- 
ations in data processing we require that at each wavelength all 
galaxies are observed by the same instrument. 

While high resolution observations are routinely available to 
resolve in detail nearby star-forming galaxies in the UV, opti- 
cal, and near-IR, the FIR is more challenging. The recent ad- 
vent of the Herschel Space Observatory (Pilbratt et al. 2010) 
now provides us with the sufficient resolution to resolve in detail 
the emission of the dust in spiral galaxies at longer wavelengths. 
The availability of Herschel data is therefore an absolute require- 
ment and fundamental to this study. The largest sample of nearby 
galaxies has been observed by Herschel in the context of the 
HRS, which observed 323 galaxies in a range of environments. 
The HRS is a K band selected, quasi-volume limited survey be- 
tween 15 Mpc and 25 Mpc. In addition, to limit the foreground 
contamination by Galactic cirrus, the sample is constituted of 
galaxies at a high Galactic latitude (\b\> 55°). 

We have selected a sample of 7 galaxies (5 in the Virgo clus- 
ter) from the 323 HRS galaxies with the following criteria: 

1 . The galaxies are normal star-forming spirals that do not have 
strong active nuclei that could strongly affect the emission in 
the UV and/or in the IR, which would induce a bias on the 
results. 

2. To perform a pixel-by-pixel analysis, the structure of the 
galaxies must be resolved at the coarsest band resolution. 
The reason is two-fold. First a good spatial resolution is 
needed to separate the various morphological components 
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Table 2. Multi wavelength data 1 . 



Bands 


Wavelength 


Instrument/Survey 


FUV 


151.6 nm 


GALEX 


NUV 


226.7 nm 




u' 


354.3 nm 


SDSS 


g' 


477.0 nm 


" 


r' 


623.1 nm 




i' 


762.5 nm 




t! 


912.4 nm 


J? 


J 


1.235 yL/m 


2MASS 


H 


1.662 yum 




Ks 


2.159 yL/m 




MIPS 70 


70 //m 


Spitzer/MIPS 


SPIRE 250 


250yL/m 




SPIRE 350 


350 yum 





(arm and interarm regions, bulge, etc.) of each galaxy and 
therefore not be only dominated by the brightest regions. 
Then, as inclination increases, different wavelengths probe 
increasingly different regions in the galaxy as UV becomes 
optically thick more rapidly than longer wavelengths, adding 
much complexity to the modelling. To constrain the FIR 
emission while preserving the highest resolution possible, 
we have decided to drop the SPIRE 500 \mv band, which has 
little consequence on the measure of the dust luminosity. To 
ensure galaxies are sufficiently resolved we have manually 
examined the images of all HRS galaxies up to 350 yum. The 
selection has been made on a case-by-case basis from this 
examination. 

3. To disentangle the emission from various components, we 
have selected only galaxies that have been observed in FUV, 
NUV, u\g\r\ i\ z\ J, H, Ks, MIPS 70 yum, SPIRE 250 yum, 
and SPIRE 350 jim. The MIPS 160 yum band was not con- 
sidered due to its coarse resolution similar to that of the 
SPIRE 500 yum band. For the FUV and NUV bands, we have 
used the GALEX (GALaxy Evolution eXplorer, Martin et al. 
2005) Nearby Galaxies Survey, Medium Imaging Survey, 
and Guest Investigator data while excluding galaxies only 
observed by the All-sky Imaging Survey as the data are 
too shallow to reach our goals. The optical data from u' to 
z' have been obtained from the SDSS (Sloan Digital Sky 
Survey, Abazajian et al. 2009). Near-IR J, H, and Ks bands 
have been acquired by 2MASS (2 Micron All Sky Survey, 
Skrutskie et al. 2006). Some of the galaxies have deeper 
near-IR images obtained by our team, but to prevent dif- 
ferences in calibration or data processing from introducing 
a bias, we have decided to use 2MASS data for all galax- 
ies. Spitzer/MIPS (Rieke et al. 2004) and Herschel/SFIRE 
(Griffin et al. 2010) data have been obtained by our team 
(Bendo et al. 2012, submitted; Ciesla et al. 2012, in prepara- 
tion) 

The detailed list of galaxies and their main parameters are pre- 
sented in Table 1. Data used from the FUV to the FIR are de- 
tailed in Table 2. 

2.2. Data processing 

To constrain the physical properties of the galaxies in a resolved 
way, a pixel-by-pixel SED from the FUV to the FIR has to be 
assembled. To do so, we have processed and convolved the data 
in a similar way as in Boquien et al. (2011). The main steps are 
described hereafter. 



Objects unrelated to the target galaxy such as foreground 
stars or background galaxies can be an important source of 
contamination when convolved to lower resolution. To limit 
this we have manually removed the brightest sources in the 
UV, the optical, and near-IR, using the imedit procedure in 

IRAF. 

For data in the UV, optical, and near-IR domains, galac- 
tic foreground extinction was corrected assuming Ry = 
3.1 with a Cardelli et al. (1989) extinction curve includ- 
ing the O'Donnell (1994) update. The differential extinction 
E(B - V) for each galaxy was obtained from NASA/IPAC 
Galactic Dust Extinction Service using the reddening maps 
of Schlegel et al. (1998). 

Processing all images to a similar point spread function 
(PSF) is crucial to study the SED in a pixel-by-pixel way 
without being affected by resolution effects. To retain as 
many resolution elements as possible while keeping strong 
constraints on the emission of the dust, we have convolved 
the data to the PSF of the SPIRE 350 fim band (24^0. We 
have used the large set of convolution kernels presented by 
Aniano et al. (2011) allowing us to degrade GALEX UV, 
optical/near-IR data as well as Spitzer MIPS data to the 
SPIRE PSE 

4. To allow for a direct pixel-by-pixel comparison, all images 
need to be projected on the same grid. For each galaxy we 
created a reference image centred on the coordinates of the 
galaxy, with a pixel size of 8'' (> 659 pc), similar to that 
of the SPIRE 350 jim image. The impact of the choice of 
the pixel size will be discussed in Sec. 4.5. This allows us to 
Nyquist sample the PSF, making it easier to distinguish the 
structures. The size of the image is taken to be twice the size 
of the circular aperture used for the integrated photometry 
presented by Ciesla et al. (2012, in preparation) to ensure 
the presence of enough sky background. All images were 
then registered on this reference image using the wregister 
procedure in iraf. 

5. For each band we subtract the background which is taken as 
the median of the pixels in an annulus with an inner radius 
of 1 . 1 times the size of the photometry aperture and with a 
width of r. The photometry aperture is taken as 1.4 times 
the optical radius as defined by Ciesla et al. (2012, in prepa- 
ration). 

6. Finally, to eliminate pixels whose physical connection to the 
galaxy is uncertain and those that are too faint, we only select 
pixels encompassed by the aperture defined by Ciesla et al. 
(2012, in preparation), and that are detected in all bands at 
a 3-cr level considering only the uncertainties on the back- 
ground level. 

An example of successive steps of the data processing from 
the FUV to the FIR is presented for NGC 4254 (M99) in Fig. 1. 

3. Spectral energy distribution fitting 

To extract the physical parameters from the observations, we 
model the SED for a large set of parameters and fit these SED 
from the UV to the FIR. We discuss the fitting procedure includ- 
ing the input and output parameters hereafter. 

3.1. CIGALE 

CIGALE is a recent code to fit the SED of galaxies simultane- 
ously from the UV to the FIR developed by our team. It has 
been used with success on nearby galaxies (Noll et al. 2009a; 
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Fig. 1. Images of NGC 4254 (M99) in FUV, NUV, g\ H, and SPIRE 250 yum (from top to bottom) at different points in data pro- 
cessing: original images (left), original images convolved to the resolution of the SPIRE 350 fim band (centre), and final convolved 
images registered to the same reference frame with a pixel size of 8'' (right). The beam size is indicated on the bottom-right hand 
of each panel. 



4 



Boquien et al.: The lRX-j3 relation on sub-galactic scales in star-forming galaxies of the Herschel Reference Survey 
Table 1. Selected sample. Data provided by Boselli et al. (2010) and Ciesla et al. (2012, in preparation). 



NGC# 


HRS# 


a 


6 


Type 


K 


Radius 


Distance 






(J2000) 


(J2000) 




(Vega mag) 


n 


(Mpc) 


3485 


33 


11 00 02.38 


+ 14 50 29.7 


Sb(r)b 


9.46 


88 


20.5 


4254 


102 


12 18 49.63 


+ 14 24 59.4 


Sa(s)c 


6.93 


258 


17.0 


4321 


122 


12 22 54.90 


+ 15 49 20.6 


SAB(s)bc; LINER; HIl 


6.59 


383 


17.0 


4519 


196 


12 33 30.25 


+08 39 17.1 


SB(rs)d 


9.56 


151 


17.0 


4535 


204 


12 34 20.31 


+08 11 51.9 


SAB(s)c;HII 


7.38 


270 


17.0 


4536 


205 


12 34 27.13 


+02 11 16.4 


SAB(rs)bc;Hll;Sbrst 


7.52 


304 


17.0 


5668 


320 


14 33 24.34 


+04 27 01.6 


SA(s)d 


11.71 


152 


22.6 



Buat et al. 2011b) and high redshift objects (Giovannoli et al. 
2011; Buat et al. 2011a; Burgarella et al. 2011) bringing new 
constraints on the physical properties of galaxies and on the at- 
tenuation laws at play. 

Briefly, CIGALE models the SED of galaxies by combin- 
ing several old and young stellar populations based on Maraston 
(2005), including nebular emission. The UV-optical domain is 
reddened by an attenuation law modified from the Calzetti et al. 
(1994, 2000) curve. This energy is re-emitted in the mid- and 
far-IR which is fitted by the Dale & Helou (2002) templates. 
The absorption and the dust emission are constrained through an 
energy balance. 

One of the main challenges in model fitting is evaluating the 
uncertainties on the main parameters. One method consists in 
finding sets of good fits by examining the ^-dimensions cube 
(in case of a standard minimisation), with n the number of 
input parameters in the model. Indeed, even if the best fit (that 
is the one that minimises the x^ value) indicates the most likely 
parameters, an arbitrarily large number of models can also pro- 
vide reasonable fits. A range of parameters can then be deter- 
mined to be acceptable. Such a method is applied in Boquien 
et al. (2010) for instance. With an increasingly larger volume 
probed by input parameters this method becomes impracticable. 
CIGALE, rather than only using a simple x^ minimisation, de- 
rives the properties of a galaxy and the associated uncertainties 
by analysing the probability distribution function for each pa- 
rameter (Walcher et al. 2008; Noll et al. 2009a). This method 
used by CIGALE to estimate parameters and their uncertainties 
is described in Kauff'mann et al. (2003). The derivation of the 
uncertainties is presented and discussed in detail in Noll et al. 
(2009a). 

3.2. Input and output parameters 

The choice of the parameters are critical. We describe here 
the parameters involved in the modelling with CIGALE, their 
range, and the accuracy of the modelling. A grid of models 
is made varying parameters determining the stellar populations 
(Sec. 3.2.1) and the conversion of UV-optical photons to the IR 
(Sec. 3.2.2). For each model of this grid the value of the x^ is 
computed and the values of the parameters are derived from the 
probability distribution function as described above. 

3.2.1. Star formation parameters 

The stellar population in galaxies is made of many generations, 
due to a complex star formation history (SFH). However, the 
relative weight of old generations is increasingly smaller in the 
UV-optical domain as they age. This creates strong degenera- 
cies in the SFH which is particularly difficult to constrain in de- 
tail. While complex SFH models give good fits, simpler models 



also give excellent fits, not allowing us to favour one model over 
the other. It is therefore reasonable to model the SED using 2 
diff'erent bursts of star formation representing the older and the 
younger stellar populations. The old stellar population is mod- 
elled with an exponentially decreasing burst over = 13 Gyr 
with an e-folding time ti ranging from 3 Gyr to 7 Gyr by steps 
of 1 Gyr, to consider the smooth, long term secular evolution 
of galaxies. On top of this burst we consider a second burst to 
model the young stars formed during the latest star formation 
episode. It has been shown on entire galaxies that a constant star 
formation rate for this episode can reproduce convincingly the 
observed SED (Buat et al. 2011b; Giovannoli et al. 2011) as the 
various, individual, short bursts of star formation are averaged 
over the entire disk of a galaxy. In the case of resolved studies, 
there are fewer individual star-forming regions averaged over in 
sub-regions in galaxies. The assumption of a constant SFR loses 
its validity as the SFH of each individual star-forming region 
becomes more important in determining the global SED of the 
galaxy sub-region. To take into account the rapid variation of the 
SFR on small spatial scales, we have chosen to use an exponen- 
tially decreasing SFH for the young stellar population. Box-like 
models have also been tested and yield nearly identical results. 
The age t2 of this burst as well as the e-folding time T2 are left as 
open parameters with 2 < t2 < 200 Myr and 1 < T2 < 100 Myr, 
each of these parameters being logarithmically spaced with a to- 
tal of 7 values each. The relative stellar mass of the young and 
the old populations, fysp, is logarithmically sampled with 6 val- 
ues in the range 0.001 < fysp < 0.316. 



3.2.2. UV-optical attenuation and infrared emission 

The presence of dust strongly aff'ects the SED in the UV and the 
optical, absorbing energetic photons to re-emit the energy in the 
IR. The distribution of the dust can strongly aff'ect the shape of 
the eff'ective dust attenuation curve (Chariot & Fall 2000; Witt 
& Gordon 2000; Panuzzo et al. 2007). In the case of a starburst, 
the slope is particularly shallow, the attenuation slowly increas- 
ing with frequency (Calzetti et al. 1994, 2000). Conversely in the 
Small and Large Magellanic clouds the extinction curve is much 
steeper (Gordon et al. 2003). In addition to intrinsic variations in 
attenuation laws, it is also well-known that young star-forming 
regions are more dusty than older ones and are therefore more 
extinguished. In the case of a starburst galaxy the young stel- 
lar population has an attenuation which is a factor ~ 2 or more 
higher than for the old stellar population (Calzetti et al. 1994; 
Wild et al. 2011). Finally, as the UV and optical emission is ab- 
sorbed and scattered, some radiation transfer eff'ects could be- 
come important. At a minimal distance of 17 Mpc, a pixel size 
of 8'' represents a physical size larger than 659 pc (1978 pc at 
the scale of the PSF). That way radiation transfer eff'ects between 
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adjacent regions should remain limited, though this strongly de- 
pends on the relative distribution of the dust and the stars. 

Changes in the attenuation laws and the presence of a UV 
bump can affect the location of the data points and the models 
in IRX-yS diagram (Burgarella et al. 2005; Boquien et al. 2009). 
In effect when only GALEX bands are available in the UV, it is 
difficult to distinguish between the presence of a bump and an 
attenuation law with a shallower slope, as the bump is contained 
in the NUV band. Indeed, in some galaxies such as the Milky 
Way the presence of a bump around 220 nm has been widely 
observed, as well as at low (Burgarella et al. 2005; Conroy et al. 
2010; Wild et al. 201 1) and high redshift (Noll et al. 2009b; Buat 
et al. 2011a). The strength of this bump strongly varies from 
galaxy to galaxy for reasons that are still subject of an intense 
debate in the literature. In addition to the effect of the metal- 
licity, recent results suggest that the inclination plays a major 
role, with edge-on galaxies exhibiting a larger bump than face- 
on galaxies (Conroy et al. 2010; Wild et al. 2011). Even though 
CIGALE can be used to detect the presence of a bump in the 
attenuation law (Buat et al. 2011a), such detailed constraints re- 
quire a fine sampling of the SED in rest-frame UV to quantify 
its amplitude accurately. In the context of the present study, the 
amplitude of the bump is expected to be weak as the selected 
galaxies are reasonably face-on (Conroy et al. 2010; Wild et al. 
2011) and fairly constant as the galaxies have similar metallic- 
ities (< 12 + logO/H >= 8.66 ± 0.11, Hughes et al. 2012, in 
preparation), close to that of the solar neighbourhood (Rudolph 
et al. 2006). While the presence of a bump mainly affects the 
NUV band, the slope of the attenuation has a broader effect 
from the FUV to the optical. To model a variation of the slope, 
CIGALE uses the starburst attenuation law as a baseline, which 
is multiplied by a factor (/l//lo)^, with A the wavelength, Aq the 
normalisation wavelength, and S the slope modifying parameter. 
Therefore, ^ = corresponds to a starburst attenuation curve, 
^ > to a shallower one, and ^ < to a steeper one, closer 
to the extinction law observed in the Magellanic Clouds for in- 
stance. To test whether the presence of a bump and a varying 
attenuation law slope have an effect on the results, we have ex- 
amined the quality of the fit with several CIGALE runs 1) adding 
a bump with an amplitude similar to that of the LMC2 supershell 
(Gordon et al. 2003) while leaving 6 free, and 2) setting the slope 
of the attenuation curve to a fixed value. Leaving the slope pa- 
rameter free brings substantial improvements to the quality of 
the fit both globally (measured with the mean;^^) and in the UV 
(measured with the mean difference between the modelled and 
the observed FUV-NUV colour). Conversely forcing the pres- 
ence of a bump slightly degrades the quality of the fits. Given 
the poor constraints, and the expectation that the bump should be 
weak, we have chosen not to consider the presence of a bump but 
we have left the slope as a free parameter ranging from S = -0.4 
to ^ = 0, with steps of 0.1, which reproduces the observations. 
We have also fixed the attenuation fraction, fan, between the 
old and young stellar populations to fau = 0.5 which is com- 
monly observed. A test run has been performed with fan = 0.75, 
showing there is little impact on the results. The attenuation of 
the young population in the V band ranges from 0.05 mag to 
1.65 mag with steps of 0.15 mag. 

The emission in the IR is handled by CIGALE using differ- 
ent sets of models or templates such as Dale & Helou (2002) or 
Chary & Elbaz (2001) for instance. For this paper we mainly 
require good estimates of the dust luminosity. As we are not 
aiming at studying the detail of the properties of the FIR emis- 
sion, the choice of templates has little consequence on the re- 
sults. Therefore we have chosen to use the Dale & Helou (2002) 



set of templates that has been extensively tested with CIGALE. 
These templates are fine-tuned to describe the properties of lo- 
cal, normal galaxies such as the ones selected in this sample. 
They are parametrised by the exponent a as defined in Dale & 
Helou (2002): dM^ (U) oc U-^'dU, with (U) the dust mass 
heated by a radiation field U. We sample this parameter between 
a = 0.5 and a = 4.0 with steps of 0.5, which covers the param- 
eter space sampled by star-forming galaxies including extreme 
cases. 

An example of a typical best-fit model found by CIGALE is 
presented in Fig. 2. 

10^ E —^—^ ■ ■ ■ I 




A[m] 



Fig. 2. Typical best fit SED by CIGALE from the FUV to the 
SPIRE 350 yum bands. This SED corresponds to one pixel in 
NGC 4254. The observed fluxes are represented by the red 
crosses with green 3-cr error bars. The value of is 0.76. The 
parameters are derived taking into account the probability distri- 
bution function yielded by a set of reasonably good fits. 



3.3. Accuracy on the output parameters 

The CIGALE code has been extensively tested by Noll et al. 
(2009a); Buat et al. (2011b); Giovannoli et al. (2011) on entire 
galaxies at low and high redshift. However, individual regions 
in galaxies are more likely to undergo local effects which can 
be explored with CIGALE. To ensure that the fitting procedure 
is reliable for individual regions in galaxies we require 1) that 
the observed SED can be reproduced by the models, and 2) that 
the intrinsic parameters are accurately estimated. To ascertain 
whether these requirements are met we apply the methods de- 
veloped in Giovannoli et al. (201 1) and Buat et al. (201 lb). 

First, to test whether the models can reproduce the observa- 
tions we compare the range of colours covered by the models to 
the colours of the observations. We concentrate on colours that 
constrain the star formation, attenuation, and properties of the IR 
emission. We compare the FUV-NUV, NUV-/, and /-MIPS 70 
colours in Fig. 3. The set of models chosen covers a larger range 
than the observations which indicates that the observations can 
be reproduced by the models. 

Another important test is to check whether CIGALE reli- 
ably estimates the output parameters. This requires an a priori 
knowledge of the intrinsic characteristics of the galaxies. One 
method is to create an artificial catalogue of galaxies, following 
Giovannoli et al. (2011). To do so we compute the best fit for 
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log F^^IF,^jy log F^.IF^^ 



Fig. 3. Comparison between the colours of the models (grey hexagons) and the observed colours (circles). For better visibility, the 
grey shade is proportional to the log of the number of models within a given bin. Black represents a high density of models whereas 
white indicates that there is no model in this bin. The colour of the circles identifies the galaxy they belong to. The left plot compares 
the NUV-r' versus the FUV-NUV colours whereas the right plots compares the FUV-MIPS 70 versus the NUV-r' colours. The 
median 3-(T error bars are displayed at the top-left corner of each plot. 



each pixel in each galaxy. This first step yields a set of artificial 
SED including their intrinsic parameters given by CIGALE. To 
take into account the uncertainties on the observations and on 
the models, for each band we add a random uncertainty drawn 
from a gaussian distribution with a taken as the observed uncer- 
tainty computed on the original images. The latter is taken as the 
uncertainty on the artificial fluxes. We then perform a new run 
of CIGALE on the artificial catalogue, using the same input pa- 
rameters as previously. In Fig. 4 we compare the parameters of 
the artificial catalogue computed from the probability distribu- 
tion function to the ones determined from the best fit along with 
the Pearson correlation coefficient p. The parameters are: 

- The stellar mass (log M^). 

- The instantaneous SFR (log SFR), the SFR averaged over 
10 Myr (log <SFR>io) and 100 Myr (log <SFR>ioo). 

- The age-sensitive Z)4000 index defined as the ratio between 
the average flux density in the 400-410 nm range and that in 
the 385-395 nm range. 

- The e-folding time of the youngest star-forming episode 

(T2). 

- The attenuation curve slope modifying parameter (6). 

- The V-band attenuation of the young (Ay^^^) and total stellar 
populations (Ay), and the FUV attenuation {Apjjy). 

- The dust luminosity (log Ldust)- 

- The a parameter as defined in Dale & Helou (2002), tracing 
the dust temperature. 

The comparison between the results from the fit and the exact 
value of each parameter shows that CIGALE yields excellent re- 
sults for most parameters with p > 0.8 for log M^, log SFR, log 
<SFR>io, Z)4000, Ay, Ay^ys p, Apuv^ and log L^ust- For 5 of them 
p > 0.9. Conversely, S is poorly constrained (p=0.363) and log 
T2 not at all (p=0.145). CIGALE yields nearly perfect estimates 
of the dust luminosity. This is due to the fact that having observa- 
tions below and beyond the peak of the dust emission constrains 
the location and the amplitude of this peak which emits the bulk 
of the dust luminosity. Not having observations at shorter wave- 
lengths has therefore only a minimal eff'ect. 



4. Results and discussion 

4.1. \RX-p diagrams and Influence of the physical 
parameters 

As mentioned in Sec. 1, the IRX-/5 diagram is one of the key 
tools to correct star-forming galaxies for the attenuation, linking 
P to IRX = log (Ldust/Lfuv)- Both terms are intimately related 
to the attenuation. 

First of all, the diff'erence between the observed UV slope 
j3 and the intrinsic UV slope in the absence of dust fio, that is 
the reddening of the UV slope by the presence of dust, can be 
directly connected to the eff'ective attenuation: 

AFuv = a^(fi-M- (2) 

This relation naturally separates the eff'ect of the SFH which en- 
tirely and exclusively determines JSq, from the eff'ect of the shape 
of the attenuation curve which entirely and exclusively deter- 
mines a^: 



(kfuv - kjsfuv) X 2.5 log Afuv/ ^nuv ' 

with k the attenuation curve. In other words, represents the 
sensitivity of the FUV attenuation to the reddening by the dust. 

At the same time IRX quantifies the relative fraction of the 
UV radiation reprocessed by dust, which permits us, under some 
approximation, to estimate the attenuation: 

A^f/y = 2.5 log (l+a/^xlO'^^), (4) 

with ajRx a constant that depends on the relative amount of emis- 
sion in the FUV band compared to attenuation sensitive bands 
(Meurer et al. 1999). Thus, combining equations 2 and 4,13 and 
IRX are linked through the following equation: 

IRX = log [[lO'-^-^^-^o) _ /^^^^j (5) 

A full derivation of equation 5 is presented in Hao et al. (2011). 

The observed slope JS is not directly accessible to us as only 
broadband UV observations are available. For a given SED, the 
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Fig. 4. Output parameters of the artificial catalogue for the best fit (x-axis) and the probability distribution function (PDF) analysis 
(y-axis) after a random error has been injected into the model SED. From the top-left comer to the bottom right corner the output 
parameters presented here are: log M^, log SFR, log <SFR>io (averaged over 10 Myr), log <SFR>ioo (averaged over 100 Myr), 
Z)4000, log T2, Ayys p, Ay, Afuv, log Ldust, and a. The Pearson correlation coefl&cient p is indicated in the top-left corner of each 
plot. For each parameter we have also computed the best linear fit which is shown in red along with the equation in the top-left 
comer where the standard deviation cr around this best fit is also indicated. These fits take into account the uncertainties computed 
by CIGALE on each data point. The median error bar computed by CIGALE is shown on the left side of each plot. 



value of /3 can vary depending on the available bands, due to 
colour eff'ects generated by the shapes of the filters. We compute 
P from the GALEX FUV and NUV bands using the following 
relation: 



log (Fnuv/Ffuv) 
log (Afuv/^nuv) 



(6) 



with F the flux density and A the wavelength. It can also be ex- 
pressed in terms of magnitudes: 



FUV -NUV 
2.5 log Ufuv/^nuv) 



■2. 



(7) 



To compute IRX = log (L^ust/Lpuv) we use the dust luminosity 
provided by CIGALE, which directly comes from the energy 
balance between dust absorption and emission. We present some 
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IRX-/5 diagrams in Fig. 5 showing the data for all the galaxies 
in order to examine the influence of various physical parameters 
estimated by CIGALE on the IRX-yS relation. 

In Fig. 5, we see that the points from the selected HRS 
galaxies lie well below the IRX-/5 starburst relations derived by 
Overzier et al. (201 1), which is in accordance to what is expected 
for non starbursting galaxies (Bell 2002; Buat et al. 2002; Kong 
et al. 2004; Gordon et al. 2004; Seibert et al. 2005; Calzetti et al. 
2005; Boissier et al. 2007; Dale et al. 2007; Johnson et al. 2007; 
Panuzzo et al. 2007; Cortese et al. 2008; Munoz-Mateos et al. 
2009; Boquien et al. 2009). There are a number of points that 
present a significantly higher IRX at a given (3 in comparison to 
the envelope described by the points. These points are roughly 
compatible with the relations defined for starburst galaxies. A 
close inspection reveals they belong to NGC 4536 which is un- 
dergoing a nuclear starburst. This shows that several diff'erent 
regimes in terms of IRX-yS can coexist within a single galaxy. 
We also see that the location in the IRX-yS diagram is strongly 
linked to some of the parameters. For instance, the stellar mass 
shows a clear gradient along the envelope. 

To understand the diff'erent trends with the parameters and 
the deviation from the IRX-yS starburst relation from the liter- 
ature, we have defined 2 quantities: the perpendicular and par- 
allel distances d^. and d\\ as shown in Fig. 6. The perpendicular 




3 

Fig. 6. Diagram showing how and d\\ are computed between 
a data point of coordinates (JS^p', IRX^p) and a given IRX-jS rela- 
tion. 

distance quantifies the deviation away from a given IRX-JS re- 
lation whereas the parallel distance quantifies gradients along 
the relation. By convention, data points that are located above 
(respectively under) the curve have d^ < (resp. > 0). To 
compute d\\ we choose the origin of the IRX-/3 curve to be set 
at IRX=-0.5. The choice of the origin has no influence on the 
results. As a reference curve we use the IRX-JS relation of Kong 
et al. (2004): IRX = log (lO^ i+o.ss^ _ ^ g^y j^e computed val- 

^ Note that a relation of the form IRX = log ( 10^^+^ - cj is eas- 
ily converted to a form similar to equation 5, with ajRx = /?o - 
(log (c) - b) /a, and = 2.5 X a. 



Table 3. Spearman correlation coeflftcients for d_i_ and d\\ versus 
various output parameters. 



Parameter 






^FUV 


0.07 


0.86 


^NUV 


0.06 


0.85 


Av 


-0.02 


0.82 


0.03 


0.76 


< Clg6 ^mass 


0.05 


0.35 


a 


-0.19 


-0.35 


6 


-0.29 


-0.62 


Z)4000 


0.04 


0.36 


XogSFR 


0.16 


0.38 


log <SFR >io 


0.15 


0.36 


\og<SFR>i^ 


0.17 


0.09 


logLdust 


0.27 


0.51 


^og Lboi 


0.33 


0.51 


logM^ 


0.29 


0.60 


logZ? 


-0.25 


-0.34 


log < b >io 


-0.24 


-0.33 


log < >ioo 


-0.11 


-0.39 



ues of the Spearman correlation coefficient of d_i_ and d\\ versus 
various output parameters are provided in Table 3 along with the 
birthrate parameter b, the ratio of the current to the average SFR 
over the lifetime of the galaxy, which is also equivalent to the 
specific SFR. The data points affected by the starbursting region 
in NGC 4536 have been discarded by only selecting points with 
IRX < 0.5/5 -h 1.5. The plots are provided in the appendix, in 
Fig. A.l and A.2 

The various parameters considered are mostly uncorrected 
with dj_. However, we see some weak structures at higher values 
of dj_ with some data points that have a higher attenuation, SFR, 
bolometric, and dust luminosities Conversely d\\ presents a con- 
vincing correlation with the majority of parameters. There is no 
trend with the birthrate parameter, either instantaneous or aver- 
aged over 10 Myr or 100 Myr. Physically, this means that none of 
these parameters seems to be directly responsible for the devia- 
tion from the starburst lRX-j3 relation, at least on a sub-galactic 
scale. 

4.2. Relation between IRX and the attenuation 

IRX and attenuation are closely linked with quantitative rela- 
tions calibrated on entire galaxies widespread in the literature 
(Kong et al. 2004; Buat et al. 2005; BurgareUa et al. 2005; 
Cortese et al. 2008; Hao et al. 2011). Whether these relations 
are also valid for individual regions in galaxies is still an open 
question. CIGALE provides us with several measures of the at- 
tenuation: Afuv, Anuv, "^Vysp^ aiid Ay. It is difl&cult to convert 
analytically from one attenuation to the other as they depend not 
only on the underlying attenuation law, but also on the relative 
luminosity of the two populations in each band. As the relation 
provided in equation 4 is not valid for longer wavelengths, fol- 
lowing Cortese et al. (2008); Buat et al. (201 lb), we have derived 
the relation between the attenuation and IRX using a fourth order 
polynomial, giving all relations the same analytic form. 



Afuv = 0.0806/ - 0.1922jc^ -h 0.6009jc^ + 0.8071x 

-hO.6621 (0- = 0.109), 
Anuv = 0.1086/ - 0.2986/ - 0.5621/ -h 0.5910jc 

-fO.4814 (0- = 0.084), 
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Fig. 5. Ratio between the IR (provided by CIGALE) and observed FUV luminosities, IRX (y-axis) versus the UV slope (x-axis). 
The colours of the individual points represent the value of the parameter indicated on the right of the colourbar. From the top-left 
corner to the bottom-right comer the parameters presented here are: log M^, log SFR, log <SFR>ioo (averaged over 100 Myr), 
Z)4000 index, Ay, Apjjy, a, S, and log <AgQ>mass- Blue points indicate a low value whereas red points indicate a high value. 
The dashed lines represent the IRX-/3 relations of Kong et al. (2004) in blue, and the Meurer et al. (1999) and Lyman-break 
analogue relations (green and red) that have been derived by Overzier et al. (2011). The median l-cr uncertainty is displayed on 
the left side of each plot. Finally, the solid black line represents the best fit for the entire sample, minus the discarded data points: 
IRX = log [(100-9W.046) _ 1^ /0.373]. 



Av^,, = -0.0382/ + 0.2322x^-0. 1732jc^ + 0.5341x 
+0.1901 (0- = 0.043), 
Ay = 0.1388/ - 0.3451/ + 0.3130/ + 0.2302x 
+0.1759 (0- = 0.046), 

where x = IRX and cr is the standard deviation around the best 
fit. Discrepant points in NGC 4536 have been discarded. 

In Fig. 7 we plot the relations aforementioned between IRX 
and the attenuation measurements. 

The relation derived for A^uv gives results close to the ones 
found in the literature for entire galaxies. Computing AAfuv as 
the diff'erence between the relation derived in this paper and the 
relation published in the literature, we find: 

- AAfuv = -0.028 + 0.032 (Burgarella et al. 2005), 

- AAfuv = 0.068 + 0.031 (Cortese et al. 2008), 

- AAfuv = 0.138 + 0.039 (Buat et al. 2011b), 

- AAfuv = 0.156 + 0.056 (Hao et al. 2011), 



The fact that we find relations close to those derived using entire 
galaxies shows that the eflfect of radiation transfer on results is at 
most minimal. 

4.3. Relation between p and the attenuation 

As mentioned earlier, assuming the unextinguished colour of the 
UV dominating population is constant, the attenuation is re- 
lated to the reddening of UV colours by equation 2: Afuv = 
a^(l3 - Po). In Fig. 8, we plot the relation between JS and the at- 
tenuation. 

We obtain the following relations assuming and JSq con- 



slant: 








Afuv 


= 0.870 OS + 2.586) 


((T = 


0.216), 


Anuv 


= 0.633 OS + 2.616) 


(CT = 


0.162), 


"^VySP 


= 0.337 OS + 2.524) 


(CT = 


0.093), 


Av 


= 0.203 OS + 2.868) 


((r = 


0.067), 
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IRX 

Fig. 7. ApuY (black points), Aj^^uv (blue points), Ay^^p (green 
points), and Ay (red points) versus IRX. The solid lines repre- 
sent the best fit with a fourth order polynomial for each attenu- 
ation measurement. The gray lines represent estimates of Apjjy 
published in the literature. Discrepant points in NGC 4636 have 
been discarded. 




-1,5 -1,0 -0.5 0.0 0.5 



Fig. 8. Apjjy (black points), Aj^jjy (blue points), Ay^^p (green 
points), and Ay (red points) versus p. The solid lines represent 
the best linear fit for each attenuation measurement after the dis- 
crepant points from NGC 4536 have been discarded. The gray 
lines represent estimates of Apjjy published in the literature. 



where cr is the standard deviation around the best fit. These 
relations do not take into account the discrepant points from 
NGC 4536. 

Compared to the relations for starburst galaxies presented in 
Overzier et al. (2011), the relations derived from the selected 
sample 1) yield a systematically lower attenuation for a given yS, 
and 2) have a shallower slope. The first point is not surprising 
as it is a well known fact that star-forming galaxies lie below 
the starburst IRX-yS relation (Kong et al. 2004), which is also 
the case for entire HRS star-forming galaxies as we will see in 
Sec. 4.6. The parameter ap in equation 2 is directly linked to the 
shape of the attenuation law (equation 3). The value = 0.870 
corresponds to a particularly extreme attenuation law. A star- 
burst attenuation law (S = 0) yields a rather grey slope of ~ 2.3, 



while S ~ -0.2 (resembling the LMC2/supershell extinction 
curve, excluding the presence of a bump) yields a slope of ~ 1.7, 
and S ~ -0.5 (resembling an SMC-like extinction curve, ex- 
cluding the presence of a bump) yields a slope of ~ 1.3. This 
shows that assuming that /3o is a constant, the shallow slope of 
the Apuy-p relation would yield an unphysically steep attenua- 
tion law. At the same time we notice that there is a considerable 
scatter around the fit. These observations suggest that might 
vary significantly across the sample depending of the actual SFH 
of each data point. We will examine this possibility in detail in 
Sec. 4.4. 

The low, unphysical value of has a large impact when 
using starburst IRX-yS relations to correct for the attenuation. 
Comparing the "M33 inner" relation obtained by Overzier et al. 
(201 1), the A/7f/y attenuation would be overestimated by 0.6 mag 
fox P = -1.5, 1.2 mag iox p - -1, 1.7 mag for - -0.5, and 
2.3 mag for P = 0, yielding errors on the SFR up to nearly an 
order of magnitude, which would overestimate the contribution 
of normal star-forming galaxies to the cosmic star formation. 

4.4. Impact of the variation ofPo 

4.4.1 . Effect of the star formation history 

In galaxies forming stars at a lower rate, the disk-averaged SFR 
can be low enough so that the intrinsic shape of the SED in the 
UV and optical bands diff'ers significantly not only from that of 
a starburst galaxy but also from object to object. Such an ef- 
fect was explored by Kong et al. (2004) in terms of the birthrate 
parameter b. Previous work (Cortese et al. 2006; Johnson et al. 
2007; Cortese et al. 2008; Boquien et al. 2009) also found a pos- 
sible eff'ect from the age of the populations. However the values 
of b, the U - B colour, the Z)4000 index, or the Ha equivalent 
width that have been used in the literature to quantify the age 
do not provide an estimate of the shape of the SED in the UV 
domain in a direct way. Indeed, they have significantly diff'erent 
timescale sensitivities and are aff'ected by attenuation at diff'er- 
ent levels. Through careful modelling of the SED, CIGALE pro- 
vides us with direct estimates of the reddening underwent by the 
intrinsic UV slope Po which allows us to correct p for the attenu- 
ation in order to retrieve Pq. In Fig. 9 (left) we present the IRX-p 
diagram with each data point colour coded according to Pq. We 
observe that there is a clear gradient of Po according to the per- 
pendicular distance from the starburst law of Kong et al. (2004) 
as confirmed in Fig. 9 (right). This strongly hints that intrinsic 
diff'erences in the UV slope play a major role to explain the de- 
viation from the starburst relation. At the same time there is also 
a non negligible scatter around the relation which suggests that 
other parameters such as the shape of the attenuation law also 
play a role. 

The shape of the attenuation curve also constrains the p- 
Afuy relation and the location of the data points in the IRX-p 
diagram (Hao et al. 2011). Indeed, a steeper slope will increase 
the reddening of the UV slope for a given quantity of dust attenu- 
ation while absorbing a larger fraction of the emission at shorter 
wavelength. Conversely, the presence of a bump in the NUV 
band will reduce the FUV-NUV diff'erential attenuation while 
only slightly increasing the fraction of reprocessed UV light into 
the FIR. As we have shown not all data points have the same Po, 
therefore a simple fit of the data points does not yield any direct 
information on the attenuation curve. At the same time, in Fig 5 
we see that data points that have a high p also tend to have a 
low S, that is a steeper attenuation curve slope. There is a visi- 
ble correlation between S and d\\ (p = -0.62). The variation of 
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Fig. 9. Left: Same as Fig. 5 with the colour of each data point corresponding to the intrinsic UV slope Po, ranging from -2.32 to 
-1.19. Data points lower than -2.32 (resp. higher than -1.19) are shown in dark blue (resp. dark red) allowing a larger dynamic 
range for the bulk of the data points. Right: Intrinsic UV slope versus the perpendicular distance d^. The Spearman correlation 
coefficient is p = 0.87. 



Table 4. Best fit Apuy = apip- in diff'erent bins of y^o- 



bin p'^ 



Po< 


-2.32 


2.28 + 0.15 


-2.28 + 0.10 


-2.32 < y 


3o<- 


2.18 


2.09 + 0.06 


-2.16 + 0.03 


-2.18 <y 


%<- 


2.04 


2.07 + 0.05 


-2.02 + 0.02 


-2.04 < / 


%<- 


1.90 


1.93 + 0.02 


-1.94 + 0.01 


-1.90 </ 


%<- 


1.76 


1.83 + 0.02 


-1.85 + 0.01 


-1.76 </ 


%<- 


1.62 


1.70 + 0.02 


-1.77 + 0.01 


-1.62 <y 


%<- 


1.47 


1.63 + 0.02 


-1.68 + 0.01 


-1.47 <y 


%<- 


1.33 


1.67 + 0.02 


-1.52 + 0.01 


-1.33 </ 


%<- 


1.19 


1.66 + 0.03 


-1.40 + 0.03 


-1.19 </?o 




1.48 + 0.05 


-1.37 + 0.05 



Pq makes it difficult to evaluate changes in the attenuation law 
across the entire sample. Therefore, to test how variations in the 
attenuation law aff'ect the estimate of Apjjy, we have divided the 
sample in 10 bins of /^o and determined the p-Afuv relation in 
each of these bins such that Afuv = [p ~ P'q)^ with ao and P'^ 
constants derived from the fit of this relation. For a fixed /^o, ffie 
relation Apjjv-P only depends on the shape of the attenuation 
curve as we saw in equation 3. The narrower range of /^o Hmits 
the eff'ects of the SFH and allows us to examine if and how y^o 
and d are linked. 

The parameters of this relation are summarised in Table 4. 
We find that data points that have Po < -1.8, that is the bluest 
UV slope, have values of ap that are similar to those of star- 
burst galaxies as determined by Overzier et al. (2011) (1.81 < 
ap < 2.07). Conversely, bins that have a higher Po tend to have 
a smaller ap, which indicates that the eff'ective attenuation law is 
steeper. The change from low to high Pq corresponds to a steep- 
ening of the attenuation curve. It can be easily understood as a 
transition between a starburst law for strongly star-forming re- 
gions that resemble starburst galaxies to attenuation laws seen 
in star-forming galaxies. As the stellar populations age, the dust 
clouds are dispersed and coherent feedback ceases. The physical 
conditions and geometry required for a starburst attenuation law 
are no longer. Finally, we notice that the value of P'^ determined 
by the fit, while close to the value of /^o obtained from CIGALE 



shows some deviations. This is likely due to variations of the 
attenuation law within each bin. 

These results confirm that the primary reason starforming 
galaxies deviate from the starburst relation is due to diff'erences 
in their intrinsic UV colour, hence a diff'erent Po. This is most 
likely due to diff'ering SFH from one galaxy to another on the 
timescale the UV is sensitive to. At the same time there is a 
clear variation of the slope of the attenuation curve, going from 
a starburst-like curve for regions in galaxies with a low Po, and 
steepening with increasing Pq. 

4.4.2. Impact of the weak constraint on 6 

As we have seen in Sec. 3.3, the constraint on d is weak. In addi- 
tion, the dynamical range on d is compressed from 0.5 to ~ 0.25 
in artificial catalogues as can be seen in Fig. 4. In turn, this com- 
pressed dynamical range exacerbates the relation between p^ and 
• To obtain strong constraints on a good sampling of the UV 
continuum with medium or narrow-band filters from the FUV to 
the U band are required (Buat et al. 2011a, Buat et al. 2012, in 
prep.), which is difficult to obtain for nearby galaxies as the UV 
is not observable from the ground. To ensure that our results are 
not aff'ected by the weak constraint on we have performed sev- 
eral test runs. First of all, we have fitted each data point setting 
6 = -0.25, which is typical of the value obtained when S is set 
free. Qualitatively we find that the results are nearly identical 
and the quality of the fits while not as good remains excellent, 
showing the major influence of Pq. Conversely, we set Po to the 
typical value of -1.7 by selecting a SFH leading to such a value. 
We also set S free to vary in a wide range well beyond realistic 
values, such as ^ < -0.5. The best fits we obtain are considerably 
worse than previously, with a significant number of them failing 
catastrophically and yielding large errors in particular in the UV 
and in the IR, sometimes over 1 magnitude. This shows that in- 
trinsically, a variation of the attenuation law is not sufficient to 
explain why normal star-forming galaxies are located under the 
starburst IRX-p curve. Conversely, variations in the SFH which 
naturally yield variations in Po can be sufficient to explain the de- 
viation. In addition, as we saw in Sec. 4.3, the relation between 
Afuv and p could only be explained by an extreme attenuation 
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law that would be even steeper than the SMC extinction law. This 
is particularly unlikely as the metallicity of the sample is close 
to the solar one. However, a spread in would naturally make 
the Apuy-p relation more shallow without requiring a variation 
of the attenuation law. 

This shows that the poor constraint on d from CIGALE does 
not really affect the results and that the variations of y^o are the 
driving reason why normal star-forming galaxies are located un- 
der the starburst IRX-/5 relation and that variations of the atten- 
uation law play only a secondary role. 

4.4.3. Relation between and the distance from the centre 

If we examine the relation between and the distance from the 
centre of each galaxy, we notice that the outer regions tend to 
have a bluer Pq than inner regions (Fig. 10). 



gradient of Po appears to be larger than what could be explained 
even by an extreme metallicity gradient. It shows that while a 
change of metallicity between inner and outer regions may play 
a role, it is not sufficient to explain the range of Pq. 

4.5. Impact of the choice of the pixel size 

The choice of a 8'' pixel size for a PSF of 24'' could influence the 
results in case of improper convolution to the lower resolution. 
To test whether our results are affected by the pixel size, we have 
reprocessed the data choosing a pixel size of 24'', and performed 
a new analysis. It turns out that the influence is weak and that our 
results are not affected. The respective envelopes in the IRX-p 
diagram and the gradients described by the two sets are similar. 
The correlation coefficients of the parameters with and d\\ 
show a variation by typically no more than 0. 1 . 
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Fig. 10. Value of Po as a function of the inclination-corrected 
distance from the centre of the galaxy in kpc. The colour of the 
circles identifies the galaxy they belong to. 

Two effects can be at play here. First there can be metallic- 
ity gradients. The nucleus of a galaxy is generally more metal- 
rich than the outer parts. The presence and the strength of these 
gradients depends on its intrinsic parameters but also whether it 
is interacting or not as radial mixing flattens metallicity gradi- 
ents (Barnes & Hernquist 1992; Kewley et al. 2010; Rupke et al. 
2010a,b). More metal-rich star-forming regions will have a red- 
der Po because of the presence of numerous absorption lines in 
the photosphere of stars. The second effect that can generate a 
bluer Po in outer regions is simply due to the presence of nu- 
merous stars in the inner regions which can contaminate the UV 
colour making Po redder. 

Disentangling these 2 effects is difficult and would require 
to have accurate metallicity maps for all the sample. To esti- 
mate the range of the effect, we have modelled a star-forming 
region with a metallicity of Z=0.02 and Z=0.001. Assuming a 
constant SFR over 10 Myr, we find that Po is ~ 0.3 dex bluer 
in the low metallicity case. The difference is even smaller when 
considering a constant SFR over a longer timescale. Such a low 
metallicity is unlikely in the large spiral galaxies in the sample, 
especially that the emission in the FIR would become particu- 
larly faint due to the depletion of dust. This gives us an upper 
bound to the expected effect. We see in Fig. 10 that for some 
of the galaxies, like NGC 4254, NGC 4321, or NGC 4535, the 



4.6. Comparison with entire galaxies 

It is unclear whether the results we have obtained on individ- 
ual data points within nearby galaxies also hold for normal star- 
forming galaxies. Indeed, if at a local level strong variations of 
Po are expected due to a quickly varying SFR, the intrinsic UV 
slope of entire galaxies is thought to undergo smaller variations 
as a much larger number of star-forming regions is averaged 
over the galactic disk compared to individual regions. This could 
for instance lead to a stronger role of a variation of the shape of 
the attenuation law. 

Being a K-band selected, volume limited sample, the HRS 
contains a large number of normal spiral galaxies. Conversely, 
IR or UV selected samples can contain a significantly larger pro- 
portion of actively star-forming galaxies that are not necessarily 
representative of normal spiral galaxies. We have selected a sub- 
set of 63 HRS galaxies that have a Sa or later morphological 
type, excluding galaxies that are too inclined requiring a/b < 3, 
with a the major axis and b the minor axis, to limit radiation 
transfer effects affecting the SED fitting. Compared to the study 
of individual regions, we forgo the SDSS u' and bands as well 
as the MIPS 70 jim one. Conversely we use ground-based U and 
V images that are deeper than their SDSS counterpart as well as 
IRAS 60 yum images as they offer a larger coverage of the HRS 
than MIPS 70 jjm images. The range of CIGALE input parame- 
ters determined for individual star-forming regions is not neces- 
sarily adapted to entire galaxies. We have chosen to use the set 
of parameters determined by Buat et al. (201 lb), with the differ- 
ence that we allow T2 to cover the same range of parameters as 
t2. That is, we allow the latest star formation episode to be over. 
Examination of the priors similar to those presented in Fig. 3 
shows that all observations of the subsample can be reproduced 
by the models. 

In Fig. 1 1 , we have plotted the IRX-p diagram, the colour of 
each galaxy representing the value of Po. We see that the range 
of IRX covered by entire galaxies is similar to the one for indi- 
vidual regions. However, some HRS galaxies have a particularly 
high p, with 19% of them having P > 0, compared to 2% for 
individual regions. The trend between Po and the perpendicu- 
lar distance that was found inside galaxies can also be retrieved 
for entire galaxies. Galaxies that have a low Po are close to the 
starburst IRX-/5 relation, and Po increases as a function of the 
perpendicular distance from the starburst IRX-p relation. We 
see that a number of galaxies also have a particularly red UV 
slope with Po > -1.2. Close inspection shows that many of these 
galaxies tend to be Hl-deficient, with the HI deficiency com- 
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Fig. 11. Same as Fig. 5. The star symbols represent the subsam- 
ple of 58 galaxies selected from the HRS that have;^^ < 5. The 
colours of individual data points and entire galaxies correspond 
to the same value of y^o, allowing for a direct comparison. 



puted as the logarithmic difference between the observed and 
the expected HI mass (Haynes & Giovanelli 1984). Their ex- 
tremely red colour is most likely due to star formation that has 
been quenched because of a lack of a gas reservoir to feed from 
(Boselli & Gavazzi 2006). We find that there is a correlation be- 
tween /^o and the HI deficiency, with p = 0.68, more gas-rich 
galaxies having a bluer /^o- 

Galaxies that have a normal HI content (HI deficiency lower 
than 0.4) have < y^o >= -1-41 ± 0.38, whereas < >= 
-0.78 + 0.72 for HI deficient ones. Such a dispersion necessar- 
ily involves important variations in the recent SFH. HI deficient 
galaxies may have lost a large fraction of their gas over a rela- 
tively short time scale due to ram pressure for instance (Boselli 
et al. 2006, 2008). Cutting abruptly star formation in such a way 
can easily create extreme values of y^o- This is in agreement with 
the result from Cortese et al. (2008) who showed that standard 
attenuation correction recipes fail for HI deficient galaxies. In 
the case of non-deficient galaxies, in addition to systematic er- 
rors, there are hints that it could be due to important variations of 
the SFH as shown by BoselH et al. (2001); Gavazzi et al. (2002); 
Boselli et al. (2009), with nearly an order of magnitude of vari- 
ation of the birthrate parameter at a given mass for a sample of 
normal star-forming galaxies. 

4. 7. The variation of the UV slope and the standard stellar 
age estimators 

The age of the stellar populations has long been suspected to be 
the reason why star-forming galaxies deviate from the standard 
starburst relations. As mentioned earlier, the birthrate parameter 
b, the Ha equivalent width, the U-B colour, or the Z)4000 index 
have routinely been used to estimate the age. However neither of 
these parameters give a good, direct estimate of JSq, but they are 
aff'ected by diff'erent specific biases. In Fig. 12 we present the 
relations between j3o, Z)4000, and the birthrate parameter. 

It turns out that neither b nor the Z)4000 index are good pre- 
dictors of JSq. The Spearman correlation coeflftcient is p = -0.52 
between logZ? and /3o, and p = 0.39 between the Z)4000 index 
and JSq. As shown in Boquien et al. (2009), this is likely due to 
mixing of successive generations of stellar populations. Indeed, 



these indicators are more accurate in case of a single, instanta- 
neous starburst which prevents contamination by other popula- 
tions. Therefore the usual age estimators tend to provide biased 
results as not all the emission is due to a recent instantaneous 
burst of star formation. The birthrate parameter is sensitive to re- 
cent star formation as well as the average star formation over the 
lifetime of the galaxy, which is closely linked to near-IR emis- 
sion but mostly unrelated to UV. Conversely, the Z)4000 index 
(or the U - B colour which straddles across the B aimer break), 
is sensitive to star formation over a longer timescale than star 
formation tracers, taking into account populations that do not 
emit significantly in the UV anymore. Simple age estimators are 
therefore not good estimators of JSq. 

5. Conclusion 

In order to understand why star-forming spirals deviate from the 
starburst lRX-j3 relation, we have modelled from the FUV to the 
FIR, including optical and near-IR data, a sample of 7 nearby, 
reasonably face-on, star-forming spirals drawn from the HRS. 
We have used the CIGALE code to estimate a large number of 
physical parameters on a pixel-by-pixel basis in each galaxy. 
The main results are: 

- The use of a canonical starburst Apuv-^ relation on normal 
star-forming galaxies may produce an overestimate of the 
SFR by almost an order of magnitude, severely hampering 
the evaluation of the contribution of these galaxies to the cos- 
mic star formation. 

- The deviation from the starburst relation cannot be explained 
in terms of the stellar mass, the bolometric or dust luminos- 
ity, the SFR (either instantaneous or averaged over 10 Myr 
or 100 Myr), Z)4000, the birthrate parameter, or the mass- 
normalised age. 

- The deviation from the starburst relation is found to be pri- 
marily due to significant variations in the intrinsic UV slope 

from one region to another. A variation of the slope 6 of 
the attenuation curve may also play a secondary role. 

- Even though the slope of the attenuation law is left uncon- 
strained by CIGALE, data points with a low value of tend 
to have a starburst-like attenuation curve whereas regions 
that have a larger p^ tend to have a steeper attenuation curve. 

- New Apjjy-p and Apuy-WJ^ relations are provided. While 
they provide statistically accurate estimates of the attenua- 
tion, due to variations of y^o and d from one region to another, 
they are not physically motivated. 

- The results obtained on a sample of star-forming spirals 
from the HRS are consistent with the ones obtained on sub- 
regions. 

These results confirm that p is not a good tracer of dust at- 
tenuation in normal star-forming galaxies if variations of the at- 
tenuation law and especially y^o are not taken into account. As 
such great care must be used when such relations are applied on 
normal star-forming galaxies without proper modelling. 
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. A.l. Parameters versus dj_. The correlation coefficient is indicated on the top-left hand corner of each plot. 
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Fig. A.2. Parameters versus d\\. The correlation coefficient is indicated on the top-left hand corner of each plot. 
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